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Abstract 

The paper continues the rigorous investigations of the mean field 
Green function solution of the effective two-dimensional two-band Hub- 
bard model [N.M. Plakida et al., Phys. Rev. B, 51, 16599 (1995)] of the 
superconducting phase transitions in cuprates, started in [Gh. Adam, 
S. Adam, J.Phys.A: Math. Theor., 40, 11205 (2007)]. Discussion of the 
{6, T) phase diagram of the model points to the divergence of the energy 
spectrum in the limit of vanishing doping 6. Finite energy spectra at all 
possible doping rates 5 are obtained provided the hopping part of the 
effective Hamiltonian is renormalized with an effective factor pointing to 
the site-pairs availability for fermion hopping processes. 

1 INTRODUCTION 

More than two decades after the discovery of the high critical tempera- 
ture superconductivity in cuprates, a theoretical model able to describe con- 
sistently the corresponding phase transition is still missing. An attractive 
candidate is the effective two-dimensional two-band Hubbard model derived 
by Plakida et al. [1] from the more cumbersome three-band p-d model [UlS] by 
use of projection techniques based on cell-cluster perturbation theories [l]-[6j 
which provide a hierarchy of the various interaction terms. 

The present paper is devoted to the scrutiny of the two-band Hubbard 
model from the point of view of its consistency with a number of essential 
features of the cuprates. The starting point of this investigation is provided 
by the rigorous results reported [7] for the generalized mean field approxima- 
tion (GMFA) solution of the thermodynamic Green function (GF) matrix of 
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the model [11 [9]. We find that, in order to get consistent solutions everywhere 
within a given cuprate family, including its reference stoichiometric structure, 
the doping rate 5 is to enter explicitly the model Hamiltonian as a renor- 
malization factor of the Hubbard 1-forms which define the fermion hopping 
conduction neighbourhoods of the spin lattice sites. 

The paper starts with the discussion of the cuprate features which are 
explicitly included into the basic model hypotheses, as well as with other 
features which have to be accommodated within the model solution to get 
a consistent description of the properties of the cuprates (section [2]). The 
Hamiltonian of the model [1] is then rewritten in an algebraically equivalent 
form (section [3|) which allows easy identification of the various contributions 
to the GMFA-GF solution. The rigorous solution of the GMFA-GF matrix 
of the model is rewritten in an alternative form which explicitly points to the 
energy spectrum divergence in the limit of vanishingly small doping rates 5 
(section . A way to secure spectrum finiteness over the cuprate families is 
discussed in section [5l The paper ends with conclusions in section [6l 

2 ESSENTIAL FEATURES OF CUPRATES 

2.1 Cu02 planes embedded in lamellar structures 

Irrespective of the specific details of their chemical composition, the 
cuprates show lamellar perovskite structures in the (a, b) plane, carrying as 
distinctive elements the Cu02 planes, separated from each other along the c- 
axis by intermediate planes or other structures of various compositions. Some 
of these intermediate structures have the role to cement the crystal structure, 
while other (charge reservoir layers in some cuprates, chain layers in other 
cuprates [lOj) provide means for the manipulation of the doping level inside 
the Cu02 planes which are generally assumed to bring the overwhelming con- 
tribution to the superconductivity in cuprates. 

Taking into account the weak connection with each other of the Cu02 
planes, the quasi-two-dimensional behaviour of the superconducting properties 
of the cuprates is simplified, within the effective two-dimensional two-band 
Hubbard model, to a single Cu02 plane. Then the effects coming from the 
weak inter-planar interaction among the Cu02 planes are incorporated via 
effective parameters into the model. 

While Geballe and Koster [TU] speculate about possible important con- 
tributions to the superconducting pairing coming from structures lying outside 
the Cu02 planes, the correct estimate of their relative magnitude cannot be 
obtained unless the quantitative Cu02 plane contributions are estimated. 
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2.2 Existence of the Fermi surface 



Fermi surface sheets in well prepared cuprate single crystals have been 
undoubtedly evidenced, first, by 2D- AGAR positron spectroscopy at the be- 
ginning of the nineties |ll|-|13j and then by ARPES and optical methods 
|14j This result is of paramount importance for the detailed description of the 
superconducting phase transition. The existence of the Fermi surface selects 
from the energy band structure contributed by the outer d-electron levels of 
the Cu ions and the outer p-electron levels of the O ions entering the Cu02 
planes only those energy bands which lay nearest to the Fermi level of the 
cuprates. 

2.3 The charge transfer insulator nature of cuprates 

The identification of the energy bands staying nearest to the Fermi sur- 
face follows from the hierarchy of the interactions resulting in energy level 
splittings and hybridizations inside these compounds. 

The crystal field effects coming from the symmetries of the crystalline 
structure at the Cu and O sites inside the Cu02 planes determine specific 
splittings of the Cu 3d levels and O 2p levels |14j . In the hole representation, 
where the vacuum state is defined by filled Cu 3d^^ and O 2p^ states, the 
resulting copper one-hole dx2_y2 state and the two oxygen one- hole Pa{x) and 
Pa{y) states inside a given elementary cell belong to a same irreducible rep- 
resentation of the point group symmetry of the cuprate. They will therefore 
hybridize among themselves resulting in characteristic energy band structures 
around the Fermi level. 

The three parameters which determine the main features of this energy 
band structure (U - the strength of the electrostatic repulsion among the 
electrons, A - the p-d interband splitting, and W ~ the band width) satisfy 
[/ > A > W , which points to the charge transfer insulator nature [TF] of the 
cuprates. Two important consequences follow. 

First, the hybridization of the copper dx2_y2 hole with the oxygen Pa{x) 
and Pain) holes results in two one- hole mixed (predominantly oxygen) p-d 
states. This interaction inside the elementary cell of the cuprate singles out 
the characteristic Zhang- Rice (ZR) singlet [16], such that the two energy bands 
lying nearest to the Fermi level are the ZR band and the upper Hubbard 
(UH) band. Thus, the effective two-band model of the superconducting phase 
transition in cuprates retains from the complete set of p-d bands inside the 
Cu02 plane precisely the ZR and UH bands and incorporates the overall small 
effects coming from the other bands into the parameters of the model. 

Second, since A ~ 2W, the resulting model Hamiltonian corresponds to 
the strong correlation limit. This leads to considerable complications in the 
mathematical characterization of the system excitations originating in inelastic 



3 



interactions since the kernel of the resulting integral representations of the 
corresponding Green functions is not separable. 

2.4 Role of the doping: cuprate families 

A high critical temperature superconductor does not simply exist as a 
well defined stoichiometric structure. It belongs to a family of cuprates (e.g., 
LSCO, YBCO, etc.), characterized by the occurrence of a characteristic stoi- 
chiometric reference structure and its specific kind of doping (either with holes, 
giving rise to a hole doped cuprate family, or with electrons, giving rise to an 
electron doped cuprate family). Within each family, the modification of the 
doping rate 5 results in drastic modifications of the physical properties. 

The reference structure is an insulator characterized by a strong anti- 
ferromagnetic exchange interaction (the highest known values of the antifer- 
romagnetic exchange parameters occur in cuprates [17J). While the antiferro- 
magnetic ordering is preserved at doping rates 6 < 6^^, a doping range exists, 
d^j. < 6 < at which a superconducting phase transition occurs with a 
doping dependent measured onset critical temperatures, Tc = Tc{6). At some 
optimum doping rate, S = 6opt G {^lrj^cr)j the critical temperature reaches 
a maximum which allows the characterization of the cuprate as a high crit- 
ical temperature superconductor. As a consequence, a given cuprate family 
exhibits a characteristic {6,T) phase diagram (see, e.g., [T ^ [T8 l [19]). 

A consistent theoretical model has to give full account of the {6, T) phase 
diagram, with correct reproduction of the antiferromagnetic ordering at zero 
and low doping rates, of the characteristic properties measured in the under- 
doped {6 < 6opt) and the overdoped {6 > 5opt) regimes respectively and, of 
course, to predict correctly the optimum doping rate 6 = 6opt at which the 
critical superconducting temperature reaches its maximum across the family. 

2.5 The effective spin lattice 

The stoichiometric reference structure {6 = 0) of a cuprate family is an 
insulator characterized by a very large gap (A ~ 2 eV [121 EO]) between the 
two subbands lying nearest to the Fermi level. As a consequence, its electron 
(hole) states are frozen at the nodes i of the two-dimensional regular lattice 
defined by the positions of the copper ions inside the Cu02 plane. 

The essentials of the behaviour of the system are preserved if the actual 
Cu02 lattice is replaced by an effective two-dimensional spin lattice having the 
spins placed at the copper sites inside the Cu02 plane. Three consequences are 
immediate for the characterization of the effective spin lattice. First, the spin 
lattice constants, ax and Uy , are given by the lattice constants of the physical 
Cu02 plane. Second, since the reference physical structure shows antiferro- 
magnetic ordering, any pair of two first order neighbouring spins within the 
reference effective spin lattice have the spins ordered in opposite directions. 
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Third, there are four such possible states at each lattice site i in the effective 
spin lattice: |0) (vacuum), \a) = {]) and \a) = ||) (single particle spin states 
inside the hole subband), and |2) = ||J,) (singlet state in the singlet subband). 

The doping of the electron states within the physical Cu02 plane is 
equivalent to the creation of defects inside the spin lattice by means of either 
spin vacancies and/or singlet states. As a consequence of the occurrence of 
spin defects, the spin lattice ceases to be frozen: hopping conduction arises. 

2.6 Hubbard operator description of hopping 

The fact that the single d electron states are tightly-bound at the sites 
of the copper ions within the Cu02 lattice accounts for the unusually low 
conduction rate in the doped cuprates. As a consequence, adequate description 
of the hopping processes between the effective lattice spin sites is got pQ in 
terms of the Hubbard operators (HOs) [21], X'^^ = \ia){if3\, where \a) and 
denote the initial, respectively final spin states at the spin lattice site i. 

At every spin lattice site i the Hubbard operator multiplication rule holds 

Xfxr = 5p,Xt\ (1) 

together with the completeness relation which secures the rigorous fulfillment 
of the Pauli exclusion principle, 

+ Xf" + Xf^ + = 1. (2) 

The single spin state creation/annihilation in a subband are described by 
fermionic HOs, while the singlet creation/annihilation, spin or charge densi- 
ties, particle numbers, are described by bosonic HOs. Therefore, the Hubbard 
operator algebra is defined both in terms of anticommutation and commuta- 
tion relations. For a pair of fermionic HOs, the anticommutator rule holds 

{Xf,x]^} = 6.,i6p,xr + Sr,o.xf), (3) 

whereas, if one or both HOs are bosonic, the commutation rule holds 

[Xf,x]^] = 6.,i6^,Xr - 6,o.xf). (4) 

Since there are two kinds of inhomogeneities (vacancies and singlets) intro- 
duced by the doping inside the spin lattice, two fundamentally different hop- 
ping conduction processes will happen: fermion hopping (of single spins to 
vacancies inside the spin lattice or single spin interband transitions) and bo- 
son hopping (of singlet spin pairs to vacancies inside the spin lattice). 

2.7 Hopping conduction neighbourhood of a spin lattice site 

Projection techniques based on cell-cluster perturbation theory |4j-[6j 
showed that the relative intensity of the hopping process relating the sites i 
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and ra of the spin lattice is determined by the non-vanishing Wannier coef- 
ficient Vim following from the overlap of the wave functions of the d-copper 
and p-oxygen states. The coefficients show (non-exponential) decrease 
with the distance rim = l^m — rj| inbetween the sites i and m. Significantly 
different from zero are the Wannier coefficients within the first three coordi- 
nation spheres around a given reference site i. An instance of typical values 
is [U [22]: for the nearest neighbouring (nn) m-sites (the first coordination 
sphere), Vim ~ z^i = 0.14; for the next nearest neighbouring (nnn) m-sites 
(the second coordination sphere), Vim ~ = — 0.13z^i, while for the m-sites 
located at the third coordination sphere, Vim ~ ^^3 = 0.16 j^i. 

The scrutiny of the hopping part of the Hamiltonian of the effective two- 
band Hubbard model ^J showed [7] (see also ^23j) that a hopping conduction 
neighbourhood of a given spin lattice site i can be defined in terms of the 
Hubbard 1-form of labels {af3,^rj), 

The actual labels (a/3, ^rj) are defined by the available in-band or inter-band 
transitions within the effective spin lattice. 

2.8 Finite hopping induced energy correction effects 

The creation of defects inside the spin lattice by doping has two kinds of 
consequences on the normal state of a cuprate. First, the emergence of hopping 
processes result in the occurrence of hopping conductivity. Second, since the 
hopping causes finite modifications of a small fraction of the electron states 
in the neighboorhoud of the Fermi surface, the correlations stemming from 
hopping induce finite corrections to the energy spectrum of the compound. 
The fulfillment of this condition by the mean field energy spectrum which 
follows from the effective two-band Hubbard model will be shown in section [5] 
to ask for a non-trivial modification of the hopping part of its Hamiltonian. 

2.9 Unconventional anomalous pairing in cuprates 

Experimental measurements resulting in inferences on the pairing mech- 
anism in cuprates concern the charge and the spin of the superconducting 
current carriers, as well as the phase of the energy gap A(k). Until the dis- 
covery of the cuprates, it was known that the anomalous pairing yielding zero 
spin and charge 2e Cooper pairs in the superconducting phase stems from 
the lattice phonon mediated interaction between conduction electron pairs in 
metallic samples. This s-wave pairing mechanism, which results in a gap func- 
tion A(k) that preserves the symmetry of the Fermi surface of the compound, 
was initially assumed in many papers to work in the cuprates as well. 
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Flux quanta measurements in cuprates at temperatures lower than Tc 
evidenced that the superconducting current is carried by electron pairs having 
the charge 2e |24j . while the small but significant drop off of the Cu Knight 
shift below Tc [25J pointed to singlet pairing in cuprates. 

Soon after the discovery of the cuprates, P.W. Anderson [26j assumed 
that the pairing mechanism should be unconventional. For the time being, 
there are several classes of models the starting hypotheses of which predict 
the occurrence of a dj.2_y2 -pairing mechanism. The two-dimensional two-band 
Hubbard model results in a static (i^2_j^2 -pairing exchange mechanism (see [9j 
and [7]). The inelastic correlation effects beyond the GMFA solution add a 
d3;2_j^2 -pairing spin fluctuation mechanism as well [^. 

This d-pairing mechanism results in a gap function A(k) the symmetry 
of which is lower than that of the Fermi surface. Phase-sensitive experiments, 
measuring the phase shifts in a dc SQUID involving a corner Josephson junc- 
tion [27] , or detecting half- flux quanta in a frustrated geometry [28] , revealed 
that the occurrence of a robust (i2.2_j^2 -pairing is a common feature of both the 
hole-doped and electron-doped cuprates. In hole-doped cuprates like YBCO, 
showing a small orthorhombic distortion from the tetragonal reference lattice 
in the Cu02 plane, a small additional s-wave component was also evidenced. 

The phase-sensitive Andreev-Saint-James spectroscopy data con- 
firmed these findings. 

2.10 Superconducting phase kinetic energy minimization 

The occurrence of the superconducting phase below some critical tem- 
perature Tc happens as a result of the minimization of the total energy of the 
system by the correlated spin configuration of Cooper pairs to a value which 
is lower as compared to that of the normal Fermi-liquid state. 

In the conventional superconductors, the electron-phonon interaction 
mechanism, which yields, below Tc, ordered Cooper pairs within a manifold of 
the electron states lying near the Fermi level, results in significant minimiza- 
tion of the potential energy of the system. 

Experimental data concerning the energy distribution of the supercon- 
ducting phase in cuprates [301 [3T] have shown that this associates the mini- 
mization of the kinetic energy of tightly bound pairs of electron states. There- 
fore, the anomalous pairing within a model consistent with these data is to 
involve essentially the kinetic energy (i.e., hopping), while vanishing or small 
correction contributions coming from the potential energy terms. 

The anomalous pairing within the two-dimensional two-band Hubbard 
model stems mainly from the hopping part of the model Hamiltonian [9, 7J, 
hence it is of kinetic origin. 
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2.11 Spin-charge separation 

The spin-charge separation in cuprates, was shown by P.W. Anderson 
|19j to provide natural understanding of the anomalous cuprate behaviour oc- 
curring in all the four distinct phases of a {6, T) phase diagram: the normal 
metallic phase, the pseudogap state separated from the previous one by the 
temperature T* = T*{6), the d-wave superconducting phase, and the anti- 
ferromagnetically ordered phase which occurs at vanishing doping or in the 
under doped regime. 

As recently discussed in f23l [33] , the spin-charge separation is recovered 
in the GMFA rigorous solution of the two-band Hubbard model ^Tj as a con- 
sequence of the fact that the spin-charge correlation matrix elements vanish 
exactly. 

3 STANDARD MODEL HAMILTONIAN 

Using ^ , the Hamiltonian of the effective two-band Hubbard model [I] 
can be rewritten in the locally manifest Hermitian form 

H = Ho + Hh, (6) 

with the single particle contribution 

i 

ho,, = EiJ2xr+E2Xf\ (7) 
and the hopping contribution 

i 

/Cii ^ / 0-0,00- Oo-,(tO\ , ^22 v-^ / 2(t,(t2 o-2,2o-\ , 
a cr 

, ^21 v-^ r, r/ 2a-,0o- 0(t,2ct\ , / crO,CT2 a2,a0\] /o^ 

(7 

In these equations the summation label i runs over the sites of an infinite two- 
dimensional lattice with the lattice constants ax and ay respectively defined 
by the crystal structure of the cuprate. The spin projection values in the sums 
over a are o" = ±1/2, a = —a. 

In 1^, El = e'd — n denotes the hole subband energy for the renormalized 
energy of a d-hole and the chemical potential fi. The energy parameter of 
the singlet subband is E2 = 2Ei + A, where A ^ Ap^ = Sp — — 2eV is an 



8 



effective Coulomb energy Ues corresponding to the difference between the two 
energy levels of the model. The hopping energy parameters fCab = '^tpdKab 
{a,b = 1,2) depend on tpd, the hopping p-d integral, and on energy band 
dependent form factors Kab- The label 1 points to the hole subband, while 2 
to the singlet subband. Inband (/Cn, /C22) and interband (/C21 = /C12) processes 
are present. 

The quasi-particle spectrum and the superconducting pairing within the 
Hamiltonian ([6]) are obtained |9j by the equation of motion technique for the 
retarded and advanced two-time 4x4 GF matrices in the (r, t)-representation, 
which represent a single matrix in the (q, u;)-representation. The GMFA so- 
lution of this matrix is summarized in the next section. 

4 MEAN-FIELD APPROXIMATION 

The Green function matrices of the model define space-time correlations 
for the four-component Nambu column operator [Sj |9] 

= (Xf2 Xp Xf (9) 

and its adjoint operator xj^ = {X'j" Xf° Xp X^''). In (P, the superscript 
T denotes the transposition. 

The retarded GF matrix is written, in Zubarev notation [32], as follows 

GiZ{t-t') = {{XUt)\xl{t'))) = -W{t-t'){{XUt),xl{t')}), (10) 

where (• • •) denotes statistical average over Gibbs grand canonical ensemble. 

The advanced GF matrix replaces in (jlOp the temporal factor by \9{t' —t). 

The GF matrix in the (r, u;)-representation is related to the GF matrix 
in the (r, t)-representation by the non-unitary Fourier transform, 

+00 

Gijait -t') = ^ J G,j^{u) e-'"(*-*') du; , (11) 
—00 

where the superscripts (r) and (a) have been omitted. 

The analytic continuations of the retarded and advanced Green func- 
tions in the complex energy w-plane define a single complex function, denoted 
Gija{'^)i with cuts (jumps) along the real energy axis. 

The energy spectrum of the Hamiltonian ([6]) is solved in the reciprocal 
space. The GF matrix in this (q, u;)-representation is related to the GF matrix 
in (r, u;)-representation by the non-unitary discrete Fourier transform 

G^3<r{^) = ^jY.^'''''^'''''^ GM,^). (12) 
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For an elemental Green function of labels (a/3, 777), we use the notation 
{{xf{t)\X]'^{t'))) in the (r, f)-representation and, similarly, {{xflX]"^))^ 
(assuming Hubbard operators at t = 0), in the (r, t(;)-representation. In the 
(q, a;)-representation, it is convenient to use the notation G"^''''''(q, w). 

We shall consider henceforth the GMFA-GF, G^{q,uj) in the form [7], 



^"(q, =x\x^- ^<x(q)l \; X= {{Xi„,xl}); 



(13) 



ia(q)=Ee'''"''^ii-; ri, =rj-r, ; Aij^ = {{[Xi„, H], xj^) . (14) 

Two kinds of particle number operators, related to the singlet subband, 
ni^ = Xr + XP, n^^ = Xr + Xf\ = n,^ + , (15) 
and to the hole subband respectively, 

nt = Xr + X^\ n% = Xr+X'l\ = nt + n% . (16) 



can be defined. The completeness relation implies njo- 



The X matrix in ()13p is diagonal 



X 



X 

6 X 



X 



X2 
XI 



= 



0' 



(17) 



where X2 and xi denote spin and site independent averages, 

X2 = M = {uia); XI = {nt) = = 1 - X2 • (18) 

In terms of the doping rate (5, it results that in the hole-doped cuprates, 

X2 = 5, XI = 1 - 5, (19) 
while in the electron-doped cuprates, 

XI = <5, X2 = 1 - <5. (20) 

To understand the consequences following for the energy spectrum, we 
rewrite the GMFA-GF p3|) in the algebraically equivalent form 



Gl{c^,uj)=x"^\l^-£M) 



x'" 



(21) 



where / denotes the 4x4 unit matrix, while <?o-(ct) is the Hermitian matrix 



4(q) = x^'^AM)x^'^ 



(22) 
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The GMFA spectrum of the model Hamiltonian ([6]) is therefore given by the 
eigenvalues of <Scr(q)- Using the results reported in [7], we get 



i?^(q) $,(q) 
($.(q))t -(^,(q)) 



4(q) = ( .r;:^.v ^t:iLt i • (23) 



The normal correlations contribute the 2x2 matrices, 

with the distinct matrix elements 

C22 = C22(q) = (^1 + A) + [022 + c^22(q)]/X2 (25) 

Cu = C7n(q) = ^1 + [022 + du{^)]/xi (26) 

C21 = C2i(q) = [a2i + d2i(q)]/VxIx^ (27) 

Here amn and (imn(q) denote respectively the one-site and two-site contribu- 
tions to the matrix w4o-(q) coming from the hopping Hamiltonian dS]): 

a22 = /Cn(rf'^°)-/C22(rf''">, (28) 

021 = (/Cn - /C22) • 2a(Tf •'^°) + /C2i((Tf ''^°) - (rf •''^)), (29) 



d 



q) = /C^nE^"^"(q)[X" + (-l)™^"x™X-] + ^'^™«x'"^q)- (30) 



mn 

0=1 



Here and in what follows, (^t^^'^"^) denotes the site-independent average of the 
Hubbard 1-form ([5]), 

{rl'n = E • ^ E(^'''^''')q7.(q) (31) 

a=l q 

for all the label sets (A//, vip) of interest. Further, 

+00 

(X^^X-^)q = ^ 1 -^i^[G^M,-^(q,u; + ie)-G^^''^^(q,u;-ie)], (32) 



while the quantities 7a (q) denote the nn (a = 1), nnn (a = 2), and third 
neighbour (a = 3) geometrical form factors, 71 (q) = 2[cos(ga;ax) + cos(gj^ay)], 
72(q) = 4cos(ga:aa:)cos(gj/ay), 73(q) = 2[cos(2ga;as) + cos{2qyay)]. 

In equation ([5n|) . {xf I = 1,2,3} denote the nn, nnn, and third co- 
ordination sphere parameters respectively coming from the phenomenological 
representation of the spin-spin correlation function (SiSj). 
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The exchange energy parameters are given by 

Jmn = 4/C„„/C2i/A, {mn} G {22, 11, 21}. (33) 
Finally, the singlet hopping contribution '^(q) is given by 

X^-'(q) = E • ^ E ^k7a(q - k) (34) 

a=l k 

where Hk = 2a{X"^ X^^)]^, while Hk = 2a{X^" X'^'')]^ for hole-doped and 
electron-doped cuprates respectively, with averages defined by (f32|) . 

The anomalous correlations contribute to (|23p the 2x2 matrices, 

' , , / -2crr2 \ , / -2aTX -Tn*! \ 

= ( -T,l 2ak ) ' ^^'^^''^^ = [ ri, 2^ J 

where 

T2 = r2(q) = [/C2261 + (1 - 5)ei?'2(q) + S^Mq)]/x2 (36) 
Ti = ri(q) = [/Cn6i + (1 - 5)ei62(q) + dCM^)]/xi (37) 
T21 = r2i(q) = [/C2161 + (1 - 5)C262(q) + <56&3(q)]/VxTxi (38) 

with the two-site exchange energies ^1 = J21 and ^2 = (<^ii + <^22)/2. 

The anomalous one-site pairing matrix elements are given by Hubbard 
1-form averages (l3T]l . 

61 = ^ 2a{Tf^^') = Y: 2a(rf '°'^) , (39) 

fj cr 

where the first expression is to be used for hole-doped cuprates, while the 
second one for electron-doped cuprates. 

The anomalous two-site pairing matrix elements following from the re- 
duction to localized Cooper pairs are 

&2(q) = E ^ E nk7«(q - k) (40) 

a=l k 

where Hk = 2a{X'''^X^'^)]^, while Hk = 2a{X^^ X^'')^, for hole-doped and 
electron-doped cuprates respectively, with averages defined in ([5^ . 

The anomalous three-site pairing matrix elements following from the 
reduction to localized Cooper pairs [7J and the splitting of the three-site terms 
as done in [33] are 

3 

&3(q)=nif)^z.,.7„(q) (41) 
a=l 
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where Ilj^''^ = 2(T(Tf ^''^^), while = 2a{T^'^'^^) for hole-doped and electron- 
doped cuprates respectively, with averages defined in (f^T]) . 

5 MODIFIED HOPPING HAMILTONIAN 

The results derived in the previous section provide the rigorous GMFA- 
GF solution for the energy matrix <?o-(q) of the effective Hamiltonian ©. 
These results have been derived assuming as starting hypotheses of the model 
six out of the eleven features discussed in section [21 namely those listed in the 
subsections ICTO and ESHO 

The features discussed in subsections 12.91 and 12.101 are immediate con- 
sequences of the results reported in Pf] and in the previous section. The 
spin-charge separation (subsection 12. lip was shown in [23ll33| to be a straight- 
forward consequence of the exact vanishing of the spin-charge correlation func- 
tions within the model. 

The features mentioned in subsections 12.41 and 12.81 ask for the finiteness 
of all the terms of the matrix i?o-(q) at any value of the doping 6, in particular 
at vanishing doping, 6 = 0. 

From the equations ([ZHD-dH]) and ([39]) -([31]), it results that both the nor- 
mal and anomalous matrix elements coming from the hopping Hamiltonian ([8]) 
are finite in the limit of vanishing doping 5 — > 0. 

Corroborating this result with the values (jl9p and (j20p of the xi and 
X2 parameters, from the equations (f25]) - (f271) and (f36]) - (f38]) it results that, in 
the hole-doped cuprates, the normal terms C22(q) and C2i(q), as well as the 
anomalous terms T2(q) and T2i(q) become infinite in the limit 5 — > due to 
the vanishing denominator X2- Similarly, in the electron-doped cuprates, the 
normal terms Cii(q) and C2i(q), as well as the anomalous terms Ti(q) and 
r2i(q) become infinite in the same limit due to the vanishing denominator xi- 

A simple remedy to this inconsistency of the standard Hamiltonian of the 
model can be proposed from the scrutiny of the reduction process resulting 
in the effective Hamiltonian ([6|). The derivation of the hopping parameters 
was done under the hypothesis of occupied d-copper and p-oxygen states of 
interest. However, under doping, part of these orbitals is empty and this fact 
is to be reflected in the occurrence of an explicit doping rate dependence of 
the effective hopping parameters. 

The simplest way is to assume the renormalization, with a convenient 
factor p, of the hopping Hamiltonian ([8]) as a whole. There are three possibil- 
ities to implement such a renormalization. 

The first is to assume p = 6. Such a hypothesis would induce, however, 
unphysical infinities in the theoretical limit 6^1. 

The second possibility is to assume p = min{xi,X2}, which would cure 
both the limits 5 — > and (5 — > 1, while inducing, however, some peculiarities 
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(turning points at 6 = 0.5) in the doping dependence of both the normal and 
anomalous matrix elements entering <?o-(q). 

The third possibility, which results in smooth dependence of the matrix 
elements of <?o-(q) on the doping rate 6 S [0, 1] is to define p = X1X2, which 
can be understood as simply assuming the site-pairs availability for fermion 
hopping processes. 

Under this hypothesis, the normal matrix elements (|25p ~ (|27|) change to 

C22 = C22(q) = (^1 + A) + [a22 + d22(q)] • Xi (42) 
Cn = Cn (q) = ^1 + ["22 + dn (q)] • X2 (43) 
C21 = C2i(q) = [021 + d2i(q)] • VXiX2, (44) 

with the singlet hopping matrix element (f3ll) replaced by 

X^-^q) = X1X2 E 4 E ^k7a(q - k). (45) 

a=l k 

The anomalous matrix elements (I36p ~ (l38p change now to 

T2 = r2(q) = [/C22fol + (1 - S)^M^) + 6CMq)] ■ XI (46) 
Ti = ri(q) = [/Cnfei + (1 - 6)^M^) + SM^)] ■ X2 (47) 
T21 = r2i(q) = [/C2161 + (1 - 5)662(q) + 5e2fe3(q)] • VXVC^, (48) 

with the isotropic one-site pairing matrix element 61, Eq. (j39p . left unchanged, 
the two-site Cooper pair term ([^U]) replaced by 

^ 1 

b2{cO = X1X2 E • ]v ^ "'^^"('1 " ^^^^ 

a=l k 

and the effective three-site term contribution (I4ip replaced by 

3 

fe3(q) = xiX2-nif)E^«-7«(q)- (50) 

It is worthwhile to note that the crystallographic symmetry of the Cu02 
lattice results in sizeable consequences on the relationships among the anoma- 
lous pairing matrix elements 61, 52(q), and 63 (q). 

In square Cu02 lattices, where = Uy (e.g. for Bi-22(n-l)n, Tl-22(n-l)n, 
Hg-12(n-l)n, and Tl-12(n-l)n cuprate families with n = 1, 2, 3), due to the 
existence of the C4 rotation axis symmetry, the anomalous averages (r{'^''^^) 
and {ti^'^'^) vanish identically. Therefore the isotropic anomalous one-site 
term 61, Eq. and the three-site term 63 (q), Eq. ([^Up . obtained after the 
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splitting |23ll33j of the three-site higher order correlation terms ([7], section 6), 
equate to zero, such that the anomalous contributions ([l6l) - (P8]) to the energy 
matrix ()23p essentially reduce to the two-site terms retained in j^. 

The rectangular Cu02 lattices (e.g., for the cuprate family YBCO), are 
very slightly different from square ones (la^ — Oyl <C mm{ax, ay}). Thus, after 
use of the C4 rotation, small non-vanishing values of the anomalous averages 
^^0-2,0-2^ and {t^'^'^'^) arise. Then the terms bi and 63 (q) in ([1H|) - (P5]) are small 
almost everywhere inside the first Brillouin zone as compared to &2(q)- Since 
the isotropic bi term points to the existence of an s-type contribution to the 
anomalous pairing, qualitative agreement exist with the experimental findings 
summarized in subsection 12.91 The three percent weight of the s-type pairing 
inferred from experiments in YBCO p3] might then be used as a constraint 
for tuning the values of the phenomenological hopping energy parameters of 
the model. 

6 CONCLUSIONS 

The scrutiny of the rigorous GMFA-GF solution of the Hamiltonian of the 
two-dimensional two-band Hubbard model of the superconducting phase tran- 
sitions in cuprates [T] unveiled the occurrence of infinite quantities, stemming 
from hopping, in the matrix elements of the energy matrix ^^(q), Eq. (j22p . in 
the limit of the vanishing doping rates 6^0. 

An analysis of the essential features of the cuprates which follow both 
from experimental data and general theoretical representations shows that the 
abovementioned infinities originate in the procedure of deriving the Hamilto- 
nian of the model, where the doping induced absence of part of the electron 
orbitals in the Cu02 plane was ignored. A phenomenological approach to the 
inclusion of this feature into the model resulted in the modification of the 
effective hopping contribution to the Hamiltonian ([6]), namely, 

H'h = XiX2Hh , (51) 

with Hh given by ([8]). Technically, this result may be viewed as a renormal- 
ization of the Wannier coefficients entering the Hubbard 1-forms dS]) with the 
factor X1X2 which expresses the availability of the spin states for fermion hop- 
ping transitions inside the hopping conduction neighbourhood of the reference 
spin lattice site i. 

This leads, in the (r, t)-representation, to expressions of the hopping 
terms of the energy matrix Sij^j which simply multiply by a factor X1X2 the 
corresponding quantities which have been obtained from Hh, Eq. ([8]). 

The rigorous reduction of the order of correlation of the boson-boson 
statistical averages involving singlet hopping (normal hopping and anomalous 
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charge-charge hopping correlations) brings supplementary X1X2 factors in the 
resulting expressions. 

As a consequence, the normal correlation terms (1230 as well as the anoma- 
lous correlation terms (j35p remain finite at any doping rates 6 both for the 
hole-doped and the electron-doped cuprates. 

The detailed investigation of the consequences of this modification of 
the Hamiltonian of the two-dimensional two-band Hubbard model will be dis- 
cussed elsewhere. 
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